Observation of Droplet Size Oscillations in a Two-Phase Fluid under Shear Flow 
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Experimental observations of droplet size sustained oscillations axe reported in a two-phase flow 
between a lamellar and a sponge phase. Under shear flow, this system presents two different steady 
states made of monodisperse multilamellar droplets, separated by a shear-thinning transition. At 
low and high shear rates, the droplet size results from a balance between surface tension and viscous 
stress whereas for intermediate shear rates, it becomes a periodic function of time. A possible 
mechanism for such kind of oscillations is discussed. 

PACS numbers: 82.70.-y, 83.60.Rs, 83.10.Gr 
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Homogenizing immiscible fluids by the use of shear flow 
is an everyday experience and a fundamental step in the 
processing of soft materials []J. By fragmenting large 
domains, the flow opposes the thermodynamics instabil- 
ity driving phase separation and leads to the formation 
of non- equilibrium steady states in which the coarsen- 
ing is stopped 0- Since the pioneering work of Taylor 
|3j on isolated Newtonian emulsions, the study of the 
rupture and deformation of isolated droplets under shear 
flow has drawn much attention. For Newtonian fluids, 
the dispersed phase forms generally somewhat deformed 
spherical droplets of the break-up size: R « x/vii where 
X is the surface tension between both phases, rj the shear 
viscosity and 7 the shear rate However, when the 
shear response of the fluids exhibits complexity as in 
semi-dilute entangled polymer solutions |4|, near criti- 
cal mixtures lamellar-sponge phase-separated mix- 
tures 0, or emulsions 0, other steady state structures 
such as cylindrical domains (known as the string phase) , 
two dimensional domains (ribbons phase) or colloidal 
crystals made of droplets are observed. Predicting the 
steady state domain morphology which results from the 
interplay between the flow, the surface tension, the vol- 
ume fraction of the two phases and their viscosities is a 
real challenge in non-equilibrium physics. Moreover, for 
a wide class of complex fluids, several different struc- 
tures can be obtained under shear flow; these steady 
states arc thus separated by out-of-cquilibrium transi- 
tions. Depending on the nature of the rheological transi- 
tion (shear-thinning or shear-thickening) and on the im- 
posed dynamic variable (stress or shear rate), structural 
bistability and/or coexistence between structures char- 
acterized by shear-banding, can be seen d, H, 0- Be- 
cause energy is constantly brought into the system, these 
transitions are also expected to lead to richer behavior 
such as bifurcations to oscillatory states or even to chaos 
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The purpose of this letter is to show for the first time 



on two-phase flow that the droplet size can be a periodic 
function of time. Those oscillations are asymptotic (i.e. 
they do not correspond to a transient regime) and dis- 
play a huge period of time, typically a few thousands of 
seconds. 

We study a pseudo binary mixture made of sodium 
bis(2-ethylhexyl) sulfo-succinate (AOT) and brine (water 
and sodium chloride) E3- At T = 25°C, for low salin- 
ities (S < 14 g.l -1 ), flat bilayers of AOT, stack upon 
each other and a lamellar phase (L a ) is observed. At 
high salinities (S > 20 g.l -1 ), the bilayers interconnect 
randomly and result in a Newtonian bicontinuous phase 
referred to in Literature as sponge phase (or L 3 phase) 
|l6j |. For intermediate salinities, coexistence between L a 
and L 3 phases is found. We prepare solutions in the two- 
phase (L a /L 3 ) region with 20% wt. AOT (from Fluka) 
and 80% wt. brine (S = 17 g.l -1 ) and let them rest for a 
few weeks to reach equilibrium. Because of a density mis- 
match between the two phases, an interface appears and 
the volume fraction of L a phase (65% at T = 25°C) can 
be measured directly from the observation of the position 
of this interface. 

Experiments: As previously described |l7| . Small An- 
gle static Light Scattering experiments (SALS) are per- 
formed at imposed shear rate with a home made trans- 
parent Couctte cell, with gap e = 1 mm and inner ra- 
dius Ri = 16 mm. In short, a circularly polarized laser 
beam (He-Ne, A = 632.8 nm) passes through the cell 
along Vw, the shear gradient direction and probes the 
sample in only one of the gaps. The scattered pattern 
corresponding to light scattered in the velocity-vorticity 
(v, z) plane is digitalized, by means of a CCD video cam- 
era for the frame acquisition. Experiments at imposed 
stress (rheology and SALS) are performed with a com- 
mercial rheomcter (RS5, Rheometrics) , using a transpar- 
ent Mooney-Coucttc cell with gap e = 0.7 mm and inner 
radius Ri = 19 mm. To investigate the effect of flow 
on the (L a /L 3 ) phase region, the solution is stirred in 
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order to obtain a macroscopic homogeneous mixture, be- 
fore introducing it into the Couette cell and shearing it 
at a constant stress. 

Results: After a transient regime which is a few hours 
long, a steady state (i.e. 7 does not any longer vary) 
is reached (see Fig. [TJ. This steady state I is charac- 
terized in SALS by a scattering ring with a well defined 
six fold modulation of its radial intensity (see Fig. GJ ■ It 
does not depend on initial conditions and is determined 
only by the value of a, the applied stress. If the flow 
is stopped, the ring persists during a few tens of min- 
utes with same size, indicating that this shear-induced 
structure is metastable. Microscopic observations reveal 
that this homogeneous structure consists of a monodis- 
perse closed-compact multilamellar droplets [13 (Inset 
Fig. |2J. The ring corresponds to the characteristic size 
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Although both steady states have similar rhcological and 
microscopic features (see inset of Fig. |5J , they can how- 
ever be easily distinguished since the scattering ring and 
therefore the characteristic droplet size vary discontinu- 
ously at the transition (Fig. |2Jl . 
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FIG. 2: Variation of R the characteristic size measured from 
the position of the scattering ring in SALS with 7. Circles 
and triangles correspond respectively to S = 17 g.F 1 and S = 
17.15 g.l - . Closed (A,») symbols represent steady states and 
open (A,o) symbols stand for the maxima and minima of the 
size during oscillations. The solid straight lines represent the 
best fits for S = 17 g.l - and S = 17.15 g.l - , having slopes: 
— 1/2. Inset: shown are the textures observed between crossed 
polarizers at 7 = 100 and 7 = 230 s -1 and the corresponding 
SALS patterns. 



FIG. 1: Steady shear rate 7 vs. stress a for T = 25° C and 
S = 17 g.l -1 (•) and S = 17.15 g.l -1 (a), the equation of 
the solid lines are (respectively for S = 17 g.l - and S = 
17.15 g.r 1 ): 7 = 0.042cr 2 and 7 = 0.036cr 2 in region I, 7 = 
0.410<T 2 and 7 = 0.620a 2 in region II. For both steady states 
(S = 17 g.l -1 and S = 17.15 g.l -1 ), n = A i7 -1/2 with Ai = 
6.05 Pa.s 1/2 and Ai = 5.18 Pa.s 1/2 , A u = 1.73 Pa.s 1/2 and 
An = L33 Pa.s 1//2 . Insets: shown are the steady droplet sizes 
(S = 17.15 g.r 1 ) for region I (■) and II (□) versus 1/cr. Solid 
lines are linear fits R oc 1/cr 

of these closed packed droplets and its modulation shows 
the existence of a hexagonal order in the positions of 
the droplets. When the stress is increased, the value of 
the steady shear rate increases accordingly to 7 oc a 2 
(indicating a shear-thinning of the sample since the vis- 
cosity 77 varies as 77 = Aj7 -1 / 2 , see Fig. QJ. The size of 
the ring becomes larger and larger indicating that the 
steady droplet size R decreases. At a well defined stress 
cr c ~ 30 Pa, a transition characterized by a jump in the 
measured shear rate occurs between steady state I and 
another steady state, labelled II. Above this transition, 
the scattering pattern consists also of a ring with six well 
defined spots and the viscosity varies as r) = Anj" 1 / 2 . 



To study the nature of the transition between states 
I and II, we now perform the same experiment at fixed 
shear rate. Figure 2 shows the evolution of the steady 
droplet size i?, deduced from SALS measurements, as 
a function of 7. Two stationary branches, labelled I 
and II are observed for respectively 7 < 100 s _1 and 
7 > 300 s . They correspond to the two steady states 
already observed in rheology (see Fig. For interme- 
diate shear rates (i.e. 100 < 7 < 300 s -1 ), very surpris- 
ingly, we do not observe any longer steady states, but 
instead oscillatory states. The size of the scattering ring 
oscillates indefinitely in time, indicating that the size of 
the droplets becomes a periodic function of time (Fig.|3J). 

Although these data have been recorded and digital- 
ized for one day (« 8 10 4 s), we have observed them visu- 
ally over more than a week, with no amplitude and period 
shifts during this time. The droplet size increases very 
slowly (the typical duration of the growth is a few thou- 
sands of seconds long) and bursts out suddenly (in a few 
tens of seconds). During the continuous growth of the 
droplets, no spatial inhomogeneities, such as bands, are 
observed in the (v, z) planes. When the laser beam prob- 
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FIG. 3: Temporal evolution of for S = 17.15 g.P 1 

and 7 = 130 s _1 . Shown is the temporal evolution of the 
scattering pattern in the (v, z) plane. 



ing the structure, is moved along the vertical z-direction, 
the size of the scattering ring (and therefore the droplets 
size) remains constant in the whole cell. During the burst 
out of the droplets, the solution becomes turbid, the 
scattering ring vanishes and horizontal bands are then 
observed in the Couettc cell. These reproducible oscil- 
lations are observed both upon increasing or decreasing 
the shear rate, for different parameters of the Couette cell 
(gaps and radii). The transition I-Oscillation occurs at 
a well defined shear rate j c m 100 s _1 whereas the tran- 
sition between Oscillation and II presents an hysteretic 
behavior for 230 < 7 < 300 s" 1 (Fig. EJ. Their period 
(a few thousands seconds) is many order of magnitude 
larger than the period of rotation of the rotor (a few sec- 
onds) and varies with 7 (see Fig.0J. These observations 
rule out any artefacts, such as a possible coupling with an 
external perturbation (building vibrations, temperature 
fluctuations...). They clearly prove that the oscillations 
correspond to a non-linear relaxation process. 

Discussion: These observations are very different from 
most unsteady flows reported in complex fluids, since 
they involve large relaxation oscillations between two dif- 
ferent structural states (I and II) 0]: tne droplet size 
oscillates. We believe that a better understanding of the 
transition I/II may provide us some insights to the os- 
cillation mechanism. For both steady states, R varies as 
j- 1 / 2 whereas rj oc 7~ 1//2 . These two scaling laws suggest 
that the droplet size results in both states, from a bal- 
ance between surface tension and viscous stress, yielding 



FIG. 4: Temporal evolution of R(t), for (a) 7 = 23 s" 1 from 
an initial shear rate of 7 = 54 s _1 , the solid line correspond 
to the best fit using an exponential expression: R = 6.08 — 
2.75exp(-7j/650), (b) 7 = 125 s"\ (c) 7 = 150 s _ \ (d) 
7 = 175 s _1 and (e) to 7 = 250 s _1 for t < 20000 s and 
7 = 280 s~ x for t > 20000 s. 

to the well-known relation Q : 

— = (1) 
777 a 

Within this picture, the surface tensions x can be ob- 
tained from the slopes of the linear fits R vs. 1/a dis- 
play in the insets of Fig. ^ This estimation respec- 
tively leads to Xi ~ 9.510 -5 N.m -1 for state I and 
X11 ~ 11.510 -5 N.m -1 for state II. Both values are com- 
patible with what is expected from a dimensional analy- 
sis, namely x ~ k B T/l 2 , where k B T is the thermal energy 
and I a correlation length of the mixture of the order 
of a few hundreds of angstroms, (typically the smectic 
distance in the phase L a or the pore size of the phase 
L 3 ). The transition I/II seems therefore related to a dis- 
continuity in the surface tension between both coexisting 
phases. We believe that a compression of the multilamel- 
lar droplets occurs at this transition yielding to release of 
the solvent outside, and as a consequence to the change 
in x- Such a phenomenon which has already been ev- 
idenced in a L a phase |20| is compatible with the dis- 
continuous jumps of the viscosity and the size observed 
at the transition I/II. It may also explain the turbidity 
enhancement occurring when the droplets burst. Within 
this picture, there is a critical pressure er c , above which 
the droplets compress and release water. More insights 
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on such mechanism can be gained by considering for in- 
stance the following differential equation, describing a 
first-order phase transition induced by Laplace pressure 
between region I (uncompressed state, (f> ~ 0) and region 
II (compressed droplets, <fi ps 1): 



i_ f xW/R 



a{2(j>- 1) + (20- 1) : 



(2) 



Next to this transition, if the shear rate is increased, the 
size R slightly decreases (since R oc 7 -1 / 2 ). Then, the 
Laplace pressure xi ' R becomes higher than o~ c and wa- 
ter is expelled from the droplets. This water release de- 
creases the viscosity and consequently the viscous stress. 
Due to Eq. the droplet size then increases again and 
the Laplace pressure decreases. As a result, the water 
goes back into the droplets and the same mechanism 
starts over leading to sustained oscillations. Because 
permeation of water through membranes is a very slow 
process [21]], such a scenario may explain the very large 
time scales involved in our observations, as suggested in 
Refs. U3. 

In a first approach, for simplicity's sake, we consider 
that the temporal evolution of R(t) is driven by the fol- 
lowing first-order kinetics: 



R = - 



1 



R 



77(^,7)7 



(3) 



with x{4>) = X/(l - 4>) + Xu<f>, and v{<t>,i) = A {4>)i 1/2 
where A((f>) = Aj(l - 0) + A u 4>. Taking for A I: A n , 
Xi, Xu-i a c the values found experimentally (Figs. ^ and 
0) and r 2 = 650 s (Fig. our simplified two equa- 
tions reproduce the two steady branches observed ex- 
perimentally for both the droplet size and the viscos- 
ity. For intermediate shear rates, in a given range of 
the a and t x parameters, it may also lead to sustained 
droplet size oscillations [22j. The transient banding be- 
havior observed during the fast burst of the droplets 
suggests that an overall description requires the set-up 
of space and time-dependent non-linear coupled differ- 
ential equations. Such approaches are currently being 
developed with success by theorists to describe spatio- 
temporal flows of complex fluids 0, 0] . We hope that 
our experiments on a model system will stimulate such 
theoretical researches. 

To conclude, we have shown in a two phase flow, the 
existence of non-linear structural oscillatory states for 
which droplet size is a periodic function of time. This 
new result raises now an important question: can the 



droplet size display chaotic behavior? We believe that 
this question will be a major topic of research in the 
study of rheochaos in the near future. 

We thank T. Douar and M. Winckert for technical as- 
sistance. We acknowledge fruitful discussions with A. 
Benayad, A. Colin, J. -P. Delville and J. Rouch. 
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